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ABSTRACT 


The probable mechanisma of fracture which may be encountered In various 
components may generally be classified in two main groups. In the first the 
fracture is of "plane strain" type which may occur in components where the 
external loads, the geometric dimensions and constraining effects, the ma- 
terial's behavior, and the environmental conditions are such that prior fo 
and during a possible fracture propagation the material is not expec :ed to 
undergo large scale plastic deformations. In this case the underlyl .g frac- 
ture theory is rather well-understood, and a criterion based on some va. Na- 
tion of (by now widely accepted concept of) fracture toughness or K . usually 
provides a highly reliable tool to deal with the problem. Thick-walled pres- 
sure vessels and other heavy-section structural components may be mentioned 
as examples which may be analyzed by using this ^articular approach. 

The second type of fracture failure which may take place in some reac- 
tor components falls into the general category of "plane stress" or "high 
energy" fracture. In a great variety of tubings and containers, due to rel- 
atively small wall thickness, large defect size, high material toughness, 
and high temperature, prior to and during a possible rupture process, around 
the defect region the material would be expected to undergo large scale 
plastic deformations. In this case the standard theories of fracture based 
on the concept of plane strain fracture toughness are not applicable. This 
type of fracture which is generally accompanied by large inelastic deforma- 
tions is (somewhat loosely) termed as the plane stress fracture for which 
currently there does not seem to be a universally accepted criterion. In 
applications to reactor components an additional complicating factor arises 
because of the fact that in this case one is dealing essentially with a 
shell of given curvature rather than a flat plate. 

The theories which are currently in use in practice to analyze plan 2 
stress type of fracture are those which are based on the concepts of crit- 
ical cracV opening stretch, K^-characterization, J-integral , and the re- 
cently pr< posed plastic instability. In this paper the application of the 
fracture criteria based on these concepts to the fracture of shells will je 
discusser and the concept of plastic instability will be developed in some 
detail. Since ch ;re is no widely accepted standard criterion to deal with 
this type of fracture which may be one of the typical failure mechanism in 
reactor componet 3 , one of the aims of the paper will be to provide an up-to- 
date critical appraisal of the current theories. 

^ ^This work was supported by NASA-Langley under the Gran' NGR-39-007-011 
and by NSF under the Grant GK-42771X. Invited paper presented at the Third 
International Conference on Structural Mechanics in leactor Technology, 

Sept. 1975, London. 
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Introduction 


In considering the fracture failure of structural solids, generally it 
has been necessary to make a distinction between two types of fractures, 
namely the so-called brittle or quasi-brittle and ductile fractures. The 
former, which usually takes place in "bulky" structural parts where the 
characteristic dimensions of the original defect causing the fracture are 
small compared to the dimensions characterizing the geometry of the part, 
is associated with relatively low fracture energy and small deformations 
prior to catastrophic failure. Ductile fracture, on the other hand, is 
associated with relatively high energy and large deformations and almost 
invariably takes place in thin plate, and shell structures. The terminolo- 
gies of flat vs. shear or plane strain vs. plane stress fracture have also 
been used to characterize the two types of fracture. Needless to say, the 
distinction is not clear cut and in practice one encounters the full spec- 
trum of failures from the highly brittle fracture of, for example, cast 
iron parts to the ductile tear of polymer sheets with varying degrees of 
inelastic deformations occuring during the process. Nevertheless, the 
loose classification has been quite useful for the purpose of identifying 
the related areas of research and developing practical fncture criteria. 

Since the spectacular examples of structural fracture failures were 
found to be the result of ductile as well as brittle type of fracture and 
partly because of the intense interest shown by the aerospace industry on 
the subject, during the early periods of the introduction of fracture 
mechanics concepts there were considerable research activities regarding 
the plane stress fracture of thin nheet structures. Emboldened by the 
success of Kj C or the fracture toughness Gj C concept characterizing the 
fracture resistance of the structural materials under plane strain condi- 
tions, at one point a single parameter characterization ot the plane stress 
fracture resistance was also thought to be possible. Thus, to determine 
the corresponding fracture resistance K^, or G^, it was thought that all one 
needs to ("o is to test a sufficiently wide panel with a large enough cen- 
tral crack. However, attempts along these lines were soor abandoned for 
the simple reason that a single parameter fracture characterization under 
"plane stress" conditions with K c as the resistance parameter did not prove 
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to bo feasible. In the Intervening years most of the efforts In the frac- 
ture research appears to have gone Into the refinement and further stand- 
ardization of plane strain fracture toughness or technology and Into 
the study of the subcritlcal fracture propagctlon phenomenon. This appar- 
ent lack of widespread Interest In research and standardize tlon attempts 
regarding ductiln fracture during this period may be due partly to the In- 
herent difficulty of the problem and partly to the fact that in most appli- 
cations usually so'iie subcritlcal crack growth under cyclic loading precedes 
the fracture process and the subcritlcal crack growth phenomenon can be 
studied very adequately by using the tools of linear fracture mechanics. 

I should be noted that a successful ductile fracture criterion re- 
quires .ot only an accurate characterization of the fracture resistance of 
the material, but also a theoretical method of evaluating a related, well- 
defined factor representing the geometry of the structure and the intensity 
of the external loads. Consider, for example, three typical ductile frac- 
ture configurations whic*’ consist of a large thin sheet with a central 
crack, a deeply edge-cracked thin strip, and a deeply edge-notched very 
thick strip under plane strain condition, all subjected to tension perpen- 
dicular to and away from the crack region. The deformation state In these 
three specimens will be entirely different. It is then intuitively clear 
that It will be difficult, perhaps even impossible to define a single factor 
which can accurately describe the Intensity of the applied load and the 
geometry of the medium at the crack tip such as, for example, the stress 
intensity factor would under conditions of plane strain fracture. The 
energy balance type of fracture criteria, however sound the underlying 
physical principles, are again a single parameter model and hence work very 
effectively only when the size and shape of the dlssip&tlcn zone around the 
crack front remain reasonably independent of the specimen geometry at the 
fracture load. During the fracture process since the (reversible and irre- 
versible) inelastic work done on the material in the dissipation zone 
absorbs a certain large percentage o£ the input energy, it is not possible 
to use .vi energy balance type fracture criterion in situations when for a 
given material the size and the shape of the dissipation zone will vary 
drastically with the geometry of the specimen, as is invitiably the case in 
materials undergoing ductile fracture. 
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Partly because of the somewhat illusive and inherently diverse nature 
of the ductile fracture process, current research regarding the material 
characterization and the development of workable fracture criteria has been 
proceeding along many different lines. In the following sections only the 
notable approaches will be discussed. Since the emphasis in the paper is 
on the review of the material from the viewpoint of applications, the dis- 
cussion will be restricted largely to the "mechanics" aspects of the problem 
and different mechanisms and models proposed r n the explanation of fracture 
growth and other materials aspects will not be considered. 

2. K^-Characterization 

2.1 The Concept and its Application 

The notion of representing the fracture resistance of thin sheet ma- 
terials by a "resistance curve"^ ^ rather than a single resistance parameter 
goes back to t..a early work done at the Naval Research Laboratory [l,2]. 

In recent years there ha3 been a considerable amount of renewed interest in 
the subject, and a great deal of research has been done on its further de- 
velopment (see the articles in [3] for a thorough review). The basis of 
the development of the concept is the observation that during the fracture 
process of thin sheet materials, depending on the specimen geometry and 
loading conditions, the unstable fracture is always preceded by a certain 
amount of stable crack extension. This is roughly due to the fact that as 
the crack length increases, because of the increasing dissipation zone size 
ahead of the crack, the resistance of the material to fracture growth also 
Increases. Thus, for a material with given thickness, as the fracture 
takes place it is possible to determine experimentally the amount of crack 
extension a-a o corresponding to a given K value. This K vs a curve (known 
as R-curve or K^-curve) may now be considered as representing the fracture 
resistance of the solid under plane stress loading conditions (for a par- 
ticular specimen geometry). Furthermore, if one can show that, or if one 
simply conjectures that the shape of this curve is independent of the in- 
itial crack length a o , the specimen geometry, and the loading conditions 

(*) 

Various other terminologies used for this purpose are: R-curve, K^-ctirve 

G^-curve, and crack extension resistance curve. 
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then one can assume that the fracture resistance of the solid under plane 
stress conditions is fully characterized by the K-curve. 


Consider, for example, the specimen and the loading condition shorn in 
Figure 1. The K^-curve shown in the figure may easily be obtained from a 
displacement controlled experiment described by Figure 2 where for a con- 
stant displacement V dK/da is always fugitive. Hence, Figure 2 represents 
a fully-stable loading configuration from which the complete K^-curve can 
be obtained. For constant load experiments shown in Figure 1 the slopes 
of loading curves K vs. a (corresponding to P»constant) are all positive. 
However, from the figure it is clear that for P<P C , it the point of inter- 
section of the loading and the resistance curves we have 


dK ^R 
da da * 


( 1 ) 


and consequently the fracture propagation will be stable. Here P c is the 
value of the load for which the loading and the resistance curves are tan- 
gent to each other. On the other hand for P«P_ and a>a , we have 

w C 


dK > dK R 


da ' da (2) 

which clearly corresponds to an unstable fri cture propagation. The critical 
crack length a and the critical stress intensity factor K„ which for the 

C l* 

given loading condition correspond to catastrophic f:i' ire are determined 
by the point of tangency of the loading (constant 1 the resistance 

(K^) curves. In an actual structure then the critical load corresponding 
to fracture instability may easily be obtained by superimposing the (calcu- 
lated) loading curves K vs. ( a-a o ) on the K^ curve and searching for the 
load level giving the tangency. For example, the O*constrnt lines shown 
in Figure 3 represent qualitatively the loading curves in a longitudinally 
stiffened panel containing a central crack. It should be emphasized that 
K-curve is assumed to be independent of a . Hence, in applications 1C- 
curve has to be translated parallel to the a axis so that its intersection 
with a axis, a^ is the same as the initial crack length in the structure. 
(For a simple graphical technique of determining K ( . sec the article by 
Creagcr In 1.3].) 

In obtaining K^-curve as well as in applications th. first question 
which has to be settled is the method of accounting for the plastic defor- 
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(nations around the crack tip in calculating the stress intensity factor K. 

This is uHual’y done by some kind of plasticity correction on the crack 

length a. In this case the alternatives are: (a) ignore the plaLtlclty 

effects and use the measured crack length a in calculations; (b) use the 

n 

so-called Irwin correction by assuming that 


a 


a < r 


ip y 




(3) 


where a is the crack length used in calculating K and 0 y is the yield 

strength; (c) use a pla-ticity correction by assuming tiiat a*a ♦ p where 

in 

p is the plastic zone size obtained from a (Dugdale-Curenblat t-type) plas- 
tic strip model; arid (d) use a compliance method to determine the adjusted 
crack length a. Method (d) is basically experimental, requires f irther in- 
strumentation during tenting (see the article k" McCabe and Heyer in [3]), 
and can be used in obtaining the K^-curve. However, it is not clear how it 
can possibly be used in calculating the loading curves K vs. a in an actual 
structure with entirely different loading and geometry. In spite of all 
its obvious limitations, largely because of its simplicity, currently 
method (b) seems to be the most widely used technique to account for plas- 
ticity correction. 


2.2 Applications to Shells 

Since the failure of thin-walled pipes, containers, and other shell 

structures containing a through crack generally falls in the category of 

plane stress fracture, it is possible to use the K_-concept to determine 

the fracture load in such structures. However, some modification of the 

current practlc" is necessary to take into account the bending and the 

curvature ef ects in shells. As a first approximation one may ignore the 

bending component of the stress Intensity factor and use the r^-plastfc- 

ity correction given by (3) wi ch K*l. , K being the membrane stress inten- 

m m 

sity factor in the shell. Tne shel, stress intensity factors are given 
usually in some numerical form (graphical or tabular) [4]. To simplify 
the applications, empirical expressions obtained through a suitable curve- 
fitting would be preferable. For example, in a cylindrical shell with an 
axial crack of length 2a, the membrane stress intensity fretur K m may be 
expressed as 
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K • A K , A • 0.481 A ♦ 0.614 ♦ 0.386e (4) 

m m p m 

where A Is Che shell parameter given by 

A - [12U-V 2 ) 14 a/^h , (5) 

R is the mean radius of curvature of the cylinder, h is the thickness, v is 
the Poisson's ratio, and <s the corresponding flat plate stress lntenlsty 
factor obtained by using che same membrane loads as in the shell (e.g., 

Kp ■ N^/rfa/h ■ p o Ri/ifa/h for a prcssurlze-J cylinder). 

For a somewhat more accurate analysis the effect of bending has to be 
taken into account and a more realistic plasticity correction must be used. 

In connection with the fatigue crack propagation in plates and shells sub- 
jected to combined bending and extension it was shown that a direct super- 
position of the stress intensity factors K m and would not produce the 
correct correlation parameter [5-7]. Thus, similar to the fatigue problem, 
one may assume that in the application of concept too the appropriate 
stress intensity factor «.ay be expressed as 

K - K ra ♦ BKh , (0<B<1) (6) 

where and are given in [4] and 6 is a constant. Selecting 0»O.5 has 
given good results in fatigue crack propagation studies, which can also be 
justified in theoretical grounds [5-7]. For the purpose of applying the 
concept one may also try the 3aice value. As for the plasticity correc- 
tion, even though one may again use the r^-correction given in (3) with K 
as expressed by (6), in this case it would be more appropriate to use the 
plasticity correction as obtained from the extension-bending strip model. 

The technique for this is described in [8] and the result.; giving the plas- 
tic zone size p maj He found in [8] and [4]. 

With regard to the acceptability of the K^-curve concept as a fracture 
criterion the main question which remains to be answered is this: are the 

basic assumptions underlying the concept, namely that the K^-curve is in- 
dependent of the initial crack length, the specimen geometry, (i.e., its 
shape and size, and the loading conditions) really valid? The experimental 
evidence regarding the results obtained from various types of specimens has 
so far been inconclusive. The comparison has been tcstrictcd almost entirely 
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Co the compact tension specimen shown In Figure 1 and the center-cracked 
tension panels. One way of 'aluatlng the effectiveness of the concept 
would be the comparison of the K f , value obtained experimentally from center- 
cracked panels with that predlctod from the K^-curve which Is found from 
the compact tension specimens. For hlgh-stiungth materials such as 7075-7& 
aluminum, T1-6A1-4V titanium alloy, and PH14-8Mo SRH1050 stainless steel 
with a simple r^-correction extremely good results have been obtained [3]. 

On the other hand the results for the high toughress materials have been at 
best "marginally unfavorable" [3]. The difficulty here of course stems 
from the fact that In the presence of very large plastic deformations the 
stress Intensity factor is no longer a realistic measure of the crack geom- 
etry and ti.e external loads and the r^-correctlon (or any plasticity correc- 
tlon> v/ould be very nearly meaningless. However, in spite of these short- 
comings l^-curve strength characterization and the related fracture cri- 
terion provide a very attractive and highly promising tool for studying 
the fracture of thin plate and shell structures containing a through crack. 

3. COD-Characterlzation 

The K^-concept described In the previous section Is applicable to thin 
plate and shell structures having through cracks only. On the other hand, 
in most thin sheet structures such as pipes, pressure vessels, and contain- 
ers the ductile fracture generally starts from defects or defect zones which 
are or may be approximated by part-through cracks. Thus, to deal with these 
as well as the through crack problems a somewhat more flexible fracture 
crlt.rlon Is needed. One such criterion is that of critical crack opening 
stretch, & cr ' The argument forming the basis of this approach is quite 
simple and is based on the assumption that in the presence of large scale 
plastic deforra.it ions the fracture process at the leading edge of the crack 
will be controlled primarily by the magnitude of local straits and the 
crack opening stretch measured or calculated at the crack fr>nt is a fairly 
good measure of these strains. Thus, according to this theo.’y the geometry 
and loading conditions and consequently the overall inelastic deformation 
state in two specimens may be quite different, but at the initiation of the 
fracture process the local conditions at the leading edge of the crack must 
have the same critical stat°. From the viewpoint of material characteriz- 
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ation the concept seems to lend itself to standardization without any diffi- 
culty [9]. In this section the theoretical results for the crack opening 
stretch in a plate of fin!-, width and in a cylindrical shell containing a 
part-through or a through crack will be presented and the correlation of a 
limited amount of data will be shown. 

3.1 Center-Cracked Plate 

Consider the plane problen described by Figure 4. It is assumed that 

a plate of width 2h and thickness b is under uniform t< nsion 0 and con- 

00 

tains r relatively large defect or a cluster of defects which may be approx- 
imated by a symmetrically located part-through crack. It will also be 
assumed that cha net ligament all around the crack is fully yielded (the 
shaded area in Figure 4b). Using the plastic strip model and replacing 
the 'trip in the yield zone by tensile tractions Oy, the problem may be 
solved by the superposition of the following three problems: 

Problem It No crack, external load: o (x,+“>) » o ; 

» yy' * O 

Problem 2: Crack: -a <;<i , v • 0: 

P P 

External load: 0 (x,0) ■ -o ,|x|<a ; 

yy o' 1 p ’ 

Problem 3: Crack: -a <x<a , y ■ 0; 

P P 

External load: o (x,0) • o v for a<|x|<* , 

yy Y ■ 1 p 

b -b 

°vv (x ’°) " "b~ °y for l x l <a * 

77 o 

Here the dlmentlons a,h,b,b (i (Figure 4 '1 and the external load O q are known 
and the length giving the plastic zone size p “ a^-a ahead of the crack 
tips is an unknown. Oy is the "flow stress" which ropres3nt3 the strain 
hardening and the yield behavior of the material and may be selected as 

Oy - (l+a)o ys , 0<a<(O u -O ys )/2o^ (7) 

where O v _ is the standard yield strength, O the ultimate strength, and ft 
an appropriately selected fixed parameter. 

Examining the loading conditions closely, it racy be seen that the 
crack problems 2 and 3 can be replaced by 
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Problem 2 t Crack: -a <x<a 

P P 


0 ; 


External load: 0 (x.0) • -o 

yy o 

• f V 1 • I'l'V 
, ° 

Problem 3 : Crack: -a <x<a , y • 0; 

P P 

External load: o^ y (x,0) ■ 0 , |x|<a , 

°yy (x,0) " b" °y • • < l x l <a p * 

This means that the part-throuRh crack problem shown in Figure 4b may be 
treated as a through crack problem provided the external load O q and the 
flow stress 0 ^ are replaced by 


o' - <3 -(1- -~-)o 


b Y 
o 


• b 

, o_ - r- a 


o b “Y 
o 


( 8 ) 


o o 

The details of the solution may be found in [ 1 0 ] . Figures 5-8 show the re- 
sults. Figure 5 gives the information to determine the plastic zone size 
a^-a. Here the parameter A^ is defined by 

A - a /h . (9) 

P P 

For Ap ■ 1 the crack plane is fully yielded and it may be shown that 

V - 1 - jr f- • (10) 

°Y b o a p 

giving the straight line in Figure 5. The curve for A»0 corresponds to in- 
finite plane for which 

TTO' 

£ 

2o’ 


a , o. 

— - cos ( — 


(ID 

P 2 °Y 

The crack opening stretch 6 calculated at the crack tip x»a (Figure 4) 
is shown in Figure 6 where the normalization factor d and the parameter A 
are defined by 


4ao Y /E , A - a/h . 


For the infinite plane A-0 and 

6 _ 2 


no 


O t , O \ 

d" -» i°* (c ° 8 ^;-> 


(12) 


(13) 


The asymp otes of the 6-curves shown in the figure correspond to the fully- 
yielded net section and are given by 
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o' 

o 




(14) 


In Che case of part-through cracks the crack opening stretch 6(x) be- 
comes maximum at x«0 which Is shown in Figure 7. For the Infinite plane 
X»0 and 6 *6(0) becomes 



^ log[(l+sinB)/cosfl] , 0 

11 



(15) 


Thus, in fracture studies based on COD considerations the relevant quantity 
representing the intensity of the external loads will be in part-through 
and 6 in through crack problems. For example, if one assumes that the frac- 
ture will start when the crack opening stretch at the crack front reaches 
a critical size 6 fr which is a characteristic strength parameter of the 
material, the load carrying capacity of the plate may be obtained from 
Figures 6 or 7 depending on whether the crack is through or part-through. 
Figure 8 shows the result for a through crack. Once 6^ r and the crack 
length a (or X«a/h) are specified the figure will give the value of O q 
corresponding to fracture initiation. 


3.2 Crack Opening Stretch Resistance Curve 

I. should be pointed out that conceived as a single parameter frac- 
ture criterion, critical crack opening stretch concept cannot accommodate 
the phenomenon of stable crack growth in thin sheet structures with a 
through crack, as pointed out in the previous section, es the crack grows 
the dissipation zone ahead of the crack and the resistance of the solid to 
fracture also grow. Consequently, to maintain the fracture propagation 
process in the plate, the stress intensity level or the rate of the exter- 
nal wc* pumped into the dissipation zone must be increased accordingly. 
Thus, particularly in the presence of large scale plastic deformations, 
since the stress intensity factor is a very poor choice to represent the 
specimen geometry and the external load, it is suggested that the crack 
extension resistance curve for sheet materials with a through crack be 
plotted by using the crack opening stretch 6 (rather thin K) as the load 
factor. In this case the experimental determination of the characteristic 
resistance curve (f^-curve) of the material and the es 1 : imtion of the 

WAN* NOV BUU2J 
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critical value of the external load or the crack opening stretch 6^ at the 
onaet of unstable fracture propagation for a given geometry and ) >" *<ng con- 
ditions would follow the same procedure as in determining K R and ^ described 
in the previous section. 


3.3 Cylindrical Shell with an Axial Crack 

Except for the effect of the shell curvature and the resulting "bulging", 
the problem for the shells is identical to the p.ate problem described in 
Section 3.1 of this paper. For a pressurized cylinder containing an axial 
crack Figure 9 shows the dimensions and orientation of the crack and the 
shell and the plastic zone size p. The details of the solution may be found 
in [8] and extensive resultB regarding the crack opening stretch are given 
in [ll]. The total crack opening stretch at any point along the crack iront 
1 ' <xpressed as 

6 t (x,z) - 6(x,0) + z0(x) , | x | <a+p , |z|< y , 

6(x,0) - v(x,+0) - v(x,-0) , 

9(x) - AiiiO). + q 2(x) # o 2 ( x) „ 2 ^ w'a.O) . (16) 

where v(x,y) and w(x,y) are, respectively, the y, and z-components of the 
displacement vector in the shell on the neutral surface, and 0(x) is the 
relative crack surface rotation (Figure 9). The following are the crack 
opening displacements of particular practical Interest: 

6 ■ 6(a,0): the conventional crack opening stretch at the crack 

3 

tips x»+a, on the neutral surface z»0. 

6 q - 6(0,0): the COD at the midpoint of the crack x-0 on the neutral 

surface. 

6 : the crack opening stretch at the midpoint x-0 and the leading 

edge z-h/2-d of a part-through external surface crack given by 

6 - 6 [1 ♦ (£ -d)/R] ♦ (5 -d)0«(O> . (17) 

c o z . < 

Figures 10-13 show some calculated sample results. Figures 10 and 11 

give 6 and 6 for a through crack. In these as well a‘» Ln the subsequent 
o a 
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figures regarding the shells the stress Oy representing the yield behavior 
of the material may again be Interpreted as a "flow stress" and is related 
to the yield and ultimate strengths of the material through (7). Other 
quantities shown in the figures are defined by 

N o - Pq R , dj - 4aO y /E , \ • [I2(l-V) 2 ]*‘«//Rh (18) 

where p q Is the internal pressure. For )»0 the problem reduces to an in- 
finite flat plate for which and 6^ are given by (15) and (13), respect- 

ively. Figures 12 and 13 show 6 o and 0 2 (O)«02 for a part-through external 
surface crack where d/h-0.5. The normalizing factor for 0^ is defined by 

d 2 - 4a0y/Eh . (19) 

If one adopts a single parameter fracture strength characterization of 
the material with the critical crack opening stretch r as the material 
coi.stant, then the load carrying capacity of the cylinder may be 'biained 
from Figures 14-18. Figure 14 shows essentially the hoop stress 
o-N /h«p R/h as a fu* tion jf the crack length as represented by X (see 

n O O 

equation 18) in a cylinder containing an axial through crack. The figure 

shows the constant 6 curves. For a given 6 and crack length the crit- 

a cr 

leal value of or p (i may be obtained by interpolation. The figure also 

give * some idea about the pressure drop necessary for crack arrest. Flg- 

u: & 5-18 show the similar results for a cylinder containing a part- 

thro' ,h external crack in which the crack opening stretch 6 at the crack 

(*) C 

front and at x-0 is the load factor of critical interest . The crack 
opening stretch 6 at the tips of a through crack of ..ame length 2a are 
shown in the figures by the small circles for d*h. Generally, the extrap- 
olated values of the hoop stress o =N /h at d*h fer constant 6 appear to 
'mailer than the hoop stress corresponding to in equel crack opening 
retch in a through crack, meaning that, at least theoretically, leak 
Defore burst l.s possible. Figure 19 shows the application of this concept 
to the results ot sjme burst testa on 2014-T6 aluminum cylinders containing 
a through crack [l3]. The figure also shows the calculated elastic stress 
intensity factor K and the plasticity corrected stress intensity factor 
Kpj, at burst pressure. K was calculated from 

^In these figures 6 C is calculated by assuming that h/R 0.465/19 [i2j. 
However, in most cases contribution of the term 6 0 z/R Ja 6 C is relatively 
small and may be neglected (see equation 17). 


( 20 ) 
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where and A^ are, respectively, membrane and bending components of the 
stress Intensity factor ratio [4,ll], p c and a c are the pressure and half 
crack length at burst. Is obtained from (20) by replacing a c by a^+p 

p being the plastic rone sire given In [8] or [4]. The figure Indicates 
that as a single fracture strength parameter the cra^k opening stretch 
appears to be preferable to the stress intensity factor with or without 
the plasticity correction. 

4. Plastic Instability 

One of the principal objections to the concept of crack opening stretch 
as a fracture criterion under conditions prevailing in ductile fracture is 
that it is a single parameter criterion and hence can only predict fracture In- 
itiation rather than unstable crack propagation. It has, for example, been 
shown that for a given material and thickness the crack initiation value 
of COD is relatively constant in a variety of specimen geometries and load- 
ing conditions (e.g., [14-15]), whereas its value at the onset of unstable 
fra ture growth may be extremely sensitive to such factors (e.g., [ 16 , 17 ]) . 
Putting aside again the microstructural factors affecting the ductile frac- 
ture process, at the continuum scale the unstable fracture may be considered 
as being the consequence of some kind of plastic instability in the liga- 
ment ahead of the crack tip. In pressurized cylinders a very simple way of 
applying this notion to obtain an estimate of the load carrying capacity 
would be trie following [ 18 3 : From the crack opening displacement results 

given by Figures 11 and 12 and similar results given in [ 1 1 3 it ma y be ob- 
served that for a through or a part-through crack, in the neighborhood of 
a certain value of the external load or the hoop stress N o /h, any small in- 
crease in the load would cause, a relatively very large Increase in 6^ or 

6 . This suggests that around this particular value of the load the phe- 
c 

nomenon taking place near the crack edge may be quite similar to the 
"necking" phenomenon observed in a ductile tensile bar wnere the material 
undergoes plastic instability. Furthermore, since the slope of the re- 
lated COD vs N /h curve increases very rapidly after a certain value of 
o 

the load, the instability load may be estimated from th< *e curves with!". 
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an acceptable degree of accuracy by assuming a fixed high slope. Figure 20 
shows the results for a cylinder with a part-through or a through crack 
assuming a slop? of approximately 20 in the normalized COD vs N^/h plots 
(Figures 11,12). The comparison of the theoretical and some experimental 
results given In [l2] on pressurized steel pipes (X52, X60V, and X60C) Is 
shown in Figures 21 and 22. In the theoretical curve shorn In Figure 21, 
as suggested in [12], the flow stress Is assumed to be Oy"O YS +10,000 psl. 

In the part-through crack case shown in Figure 22 the experimental results 
were obtained for 0.492<d/h<0.5ll and the theoretical curve is based on 
d/h*0.5 and O Y *(l-K).Ob)0 Y g. It may be observed that the agreement seen In 
Figures 21 and 22 appears to be rather good. However, it should also be 
pointed out that because of the presence of slow stable crack growth pre- 
ceding the onset of unstable fracture, the correlation of these experimental 
results with a critical crack opening stretch as the fracture criterion was 
not as successful. 

The results given in [12] indicate that for near-failure conditions in 
pressurized cylinders with a part-through crack generally the era .k opening 
displacement 5 at the mid-section x«0 (see equation 17) is greater than 
O.l(h-d). On the other hand, investigations of the edge-cracked specimens 
by the J-integral method [19] suggest that for the appl lrability of this 
(J-integral) method of characterization, 6 c needs to be less than 0.04(h-d). 
Thud, it is clear that relative to 6 c tie thicktess of the net ligament is 
not large enough to characterize its failure in terms of a progressive 
crack extension model and, since the net ligament is fully-yielded, some 
alternative model based on plastic instability is needed. Such a model is 
described in [ L8] . 

In the case of a "deep" part-through axial crack if the pressure is 
such that the net ligament is fully yielded, then within the net ligament 
the strains may approximately be expressed as 

e-0,e+e-0 (21> 

* y z 

Thus, a nominal value of the average tensile strain in the net ligament may 
be estimated by 

e - 6 e /(h-d) (22) 
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where 6 ( is the average net ligament stretch evaluated at the midpoint of 
the net ligament and is given by (see (16)) 

6, - 6(0, -d/2) - 6 (l-d/2R) - 0-(0)d/2 . (23) 

C 0 2 

For the steel pipe experiments reported in [12] near failure conditions e 
is estimated to be above 0.15 which is slightly less than the expected max- 
imum strain in an unnotched tensile bar but greater than that for a deeply 
notched bar. To illustrate the stability behavior of the net ligament a 
model can be construe*-®' 4 the following manner. Let us assume that for 

any € greater than, say, 0.1 the effective crack depth may be expressed as 

d ■ d + a(e)6 (24) 

o e 

where d Q is the Initial crack depth (prior to loading) and a (c) is a co- 
efficient (0<a<l) approaching unity as e increases into the range of 0.3 to 
0.5. If we assume that a(c) is a known function and note that for given R, 
h,a, and p 6 is a function of d, then, with (22), (24) provides a (highly 

O b 

nonlinear) equation to determine the value of d. This can be done by a 
successive approximation scheme as follows: 

d(0) - d , d(N) - d + a(e(N-l))6 (N-l) , N - 1,2,... (25) 

O O E 

It is clear that if this successive approximation converges, then as 
a result of the load redistribution there will be a stable crack configura- 
tion. On the other hand the divergence of the Iteration in (25) will imply 
a net ligament plastic instability. The basic mechanics of this stability 
model will depend, among other factors, on the selection of a(e) and the 
flow stress o y . Using constant values of a between 0.4 nnd 0.8 it was round 
that roughly equivalent failure conditions can be predicted with small ot 
and o y well above o ys or large a and a y close to o ys [lb]. 

In the Battelle experiments described in [ 1 2 ] the outer surface COD 
defined by 

6 o - 6(0, h/2) » 6 q (1 + | 6 2 (0) (26) 

vas also measured. The (average) measurements and other relevant informa- 
tion were: 

R - 18 in., h - 0.403 in., d Q = 0.201 in., 2a « 3.8 in., o ys - 64.6 ksi 
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P o (ksi): 1.0 1.2 1.25 1.29 

6 (ailU)s 29 60 80 failure 

The succ live approximation scheme va« used to evaluate 6^/h by assuming 
that [18] 

°Y “ °YS + 2 k8t » a(e) * °‘ 46 + 0.54(1 ~) . (27) 

The results are shown in Figure 23. In the figure the curves, A,B,C, and D 
correspond to the pressures 1.0, 1.2, 1.25, and 1.267 ksl, respectively. 

It Is seen that 6 q for the first three pressures converges and for p q ■ 

1.267 ksl It diverges. Increasing with N In nearly linear fashion and corre- 
sponding to a net ligament Instability. The figure also shows the calcu- 
lated values of 6^ at stability (and for N“14 for curve D) . Even though 
the iteration scheme predicts the failure pressure qultr accurately, it 
predicts Increasingly smaller values for 6 at lower pressures. These re- 
sults seem to be encouraging enough to warrant further investigation of the 
plastic instability concept along the lines described in this section. 

5. A Two-Parameter Fracture Criterion 

As indicated before, a single parameter fracture criterion is inade- 
quate to characterize the ductile frarture particularly in thin sheet 
structures. The resistance curve concept based on any of the load factors 
such as G, K, 6, or J in this respect may be considered, at least theo- 
retically, as having infinitely many discrete parameters in the material 
characterization. From the practical viewpoint, the shortcomings of the 
concept are obvious, namely that it is purely empirical, requires exten- 
sive experimental work for characterization, and is not easy to apply. It 
is therefore natural to think in terms of a multi-parameter criterion as a 
practical alternative. Such a modil based on two parameters has been de- 
veloped by Newman [ 20 , 21 1 . 

The idea of the two parameter model is partly based on the nonlinear 
or elastic-plastic stress and strain concentration at the root of a blunt 
notch or a crack considered by Neuber [22] and others [23-27]. Fi.'it, 
Neuber's result is expressed in the form [20, 2l] 

oeE “ o 2 (28) 

e 
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whe *e O and e are the local stress a..d strain at the notch root, E Is the 
"modulus" and 0 Is the local elastic stress. Essentially, the asymptotic 
expression for the crack tip stress with two terms In the expansion Is then 
substituted into (28), leading to, after some manipulations, an expression 
of the form 


le 


"3n 


« •• 

i - m — 

*u 


(S <S ) 
n u' 


(29) 


where and m are the two material parameters (characterizing the fracture 
strength), K le is the linear elastic stress Intensity factor, S n is the 
applied stress and is the maximum possible nominal elastic stress com- 

puted from the load which would produce the ultimate strength 0^ on the 


K. Is calculated from the linear elastic fracture 
le 


complete net section, 
mechanics. 

dltions and are some multiples of the magnitude of the external load and 


S n and S u depend on the specimen geometry and the loading con- 


the ultimate strength O u , respectively. 


For a given material of given 


thickness and m are determined from the experimental fracture test re- 
sults by using (29) and a least square method [20]. Note that in (29) the 
denominator represents the effect of material nonlinearity. For m*0 the 
equation reduces to the linear elastic fracture model used for plane strain 
problems. For m»l the model Is similar to that described in [28] for high 
toughness materials. 


The model given by (29) has been used in [20] and [21 ] to study and 
correlate the results of fracture tests on thin sheets of a relatively 
good variety of materials ai.u specimen geometries having either a part- 
through or a through crack. The tesult appeavs to be extremely successful 
Indicating that there is a great deal of promise in trying to characterize 
the ductile fracture by a multi-parameter model. 


6. J-Integral 

The J-integral which was originally developed for nonlinear elastic 
(i.e., non-disslpat ive) materials [29] has now become, along with K and 6, 
one of t;;e standard load factors (representing the magnitude of the 
applied load and the specimen geometry) In fracture studies. As a basis 
of a ductile fracture criterion thin concept too suffers from all the 
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disadvantages of a single-parameter material characterization model. It 
is further restricted to two-dimensional crack geometries and as yet cannot 
be used In such problems as shells, part-through cracks, and other three- 
dimensional crack geometries. However, together with the concept of a re- 
sistance curve or a multi-parameter strength characterization it may prove 
to be an effective tool in ductile fracture studies for thin sheet struc- 
tures having a through crack. 

7. Conclusions 

(a) From the viewpoint of structural failure^ ductile fracture process 
has two aspects, namely the fracture initiation followed by a stable crack 
growth and the onset of unstable fracture propagation. Fracture initiation 
occurs when the local conditions at the leading edge of the crack reach a 
critical state and hence, in principle, can be predicted by a suitable 
single parameter criterion such as crack opening stretch or, in plane crack 
geometries, J-tategral. 

(b) The onset of unstable crack propagation is basically a stability 
phenomenon which is controlled by both the local conditions at the leading 
edge of the crack and the global state of energy flow Into the dissipation 
zone ahead of the crack front. Because of the large size of the dissipa- 
tion zone and the plastic strains, the energy ilow process in ductile 
fracture propagation is a highly complex phenomenon. As the crack propa- 
gates some of the (externally added or internally released) input energy is 
transformed into the energy of residual stresses stored In the wake left 
behind the crack, there is certain amount of dissipation due to heating 
and restructuring the material, and there may be some form of surface 
energy. Nevertheless, the fracture instability is almost certainly tre 
result of the rate of input energy exceeding that of stored and dissipated 
energies, with the critical local fracture conditions aj neccesary condi- 
tions. This being the case, clearly a single paraneter fracture criterion 
is not adequate to characterize the phenomenon. 

(c) The ductile fracture propagation process may be characterized bv 
either a multi-parameter (discrete) model or some type of a "resistance 
curve" which may be considered as a continuous model expressed graphically. 
Some of the existing models have been described in this niper. These 
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models may have certain shortcomings but they do Indicate the direction of 
necessary research In the field. 

(d) There Is a group of problems (notably the part-through crack prob- 
lems in plates and shells) In which the ductile fracture process cannot be 
modeled as a progressive crack growth pnenomenon. In such cases a "net 
ligament plastic Instability" type of model appears to be adequate to char- 
acterize the fracture process. 
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Figure 5. Master curves for determination of plastic rone size in 
cracked strip with fully-yielded net ligament. 


• / 




Figure 7. Crack opening stretch at the center x«0 (d = 4aOy/E) . 






Figure 10. Crack opening displacement 6 in the neutral surface z*0 and 

_ J * ^ J /. _ C -1 0_l 1 1 I A.L. J _ 1 A. 1. 


at the midsection x»0 of the shell with an axial through crack 
<d x - 4ao Y E , X - [12(1^2)]^ a/v^h) . 
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esc N /h vs. A for a constant crack 


,0) in a cylindrical shell with an 
o y / E , A - [12(1-V Z )P a/®. 
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Figure 21. Comparison of the experimental [12] and theoretical results for 
pressurized steel cylinders with an axial throug?. crack. 




Figure 23. Successive approximation values of 6 /h for the Battelle Expert 
ment No. 70-13 [12,18]. <N is the number of iterations.) 




